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ABSTRACT 

Activity of the nuclei of galaxies and stellar mass systems involving disk accretion to black holes is 
thought to be due to (1) a small-scale turbulent magnetic field in the disk (due to the magneto-rotational 
instability or MRI) which gives a large viscosity enhancing accretion, and (2) a large-scale magnetic 
field which gives rise to matter outflows and/or electromagnetic jets from the disk which also enhances 
accretion. An important problem with this picture is that the enhanced viscosity is accompanied by an 
enhanced magnetic diffusivity which acts to prevent the build up of a significant large-scale field. Recent 
work has pointed out that the disk's surface layers are non-turbulent and thus highly conducting (or non- 
diffusive) because the MRI is suppressed high in the disk where the magnetic and radiation pressures are 
larger than the thermal pressure. Here, we calculate the vertical (z) profiles of the stationary accretion 
flows (with radial and azimuthal components), and the profiles of the large-scale, magnetic field taking 
into account the turbulent viscosity and diffusivity due to the MRI and the fact that the turbulence 
vanishes at the surface of the disk. We derive a sixth-order differential equation for the radial flow 
velocity Vr{z) which depends mainly on the midplane thermal to magnetic pressure ratio /3 > 1 and the 
Prandtl number of the turbulence V = viscosity/diffusivity. Boundary conditions at the disk surface 
take into account a possible magnetic wind or jet and allow for a surface current in the highly conducting 
surface layer. The stationary solutions we find indicate that a weak (/3 > 1) large-scale field does not 
diffuse away as suggested by earlier work. For a wide range of parameters /3 > 1 and P > 1, we find 
stationary channel type flows where the flow is radially outward near the midplane of the disk and radially 
inward in the top and bottom parts of the disk. Channel flows with inward flow near the midplane and 
outflow in the top and bottom parts of the disk are also found. We find that Prandtl numbers larger 
than a critical value (estimated to be 2.7) are needed in order for there to be magnetocentrifugal outflows 
from the disk's surface. For smaller V, electromagnetic outflows are predicted. 

Subject headings: accretion, accretion disks 
binaries 



galaxies: jets — magnetic fields — MHD — X-rays: 



1. INTRODUCTION 



Early work on disk accretion to a black hole argued that 
a large-scale poloidal magnetic field originating from say 
the interstellar medium, would be dragged inward and 
greatly compressed near the black hole by the accret- 
ing plasma (Bisnovatyi-Kogan & Ruzmaikin 1974, 1976) 
and that this would be important for the formation of 
jets (Lovelace 1976). Later, the importance of a weak 
small-scale magnetic field within the disk was recognized 
as the source of the turbulent viscosity of disk owing to 
the magneto-rotational instability (MRI; Balbus & Haw- 
ley 1991, 1998). Analysis of the diffusion and advection of 
a large-scale field in a disk with a turbulent viscosity com- 
parable to the turbulent magnetic diffusivity (as suggested 
by MRI simulations) indicated that a weak large-scale 
field would diffuse outward rapidly (van Ballegooijen 1989; 
Lubow, Papaloizou, & Pringle 1994; Lovelace, Romanova, 
& Newman 1994, 1997). This has led to the suggestion 
that special conditions (nonaxisymmetry) are required for 
the field to be advected inward (Spruit & Uzdensky 2005). 
Recently, Bisnovatyi-Kogan and Lovelace (2007) pointed 
out that the disk's surface layers are highly conducting 
(or non-diffusive) because the MRI is suppressed in this 
region where the magnetic energy-density is larger than 
the thermal energy-density. Rothstein and Lovelace (2008) 



analyzed this problem in further detail and discussed the 
connections with global and shearing box magnetohydro- 
dynamic (MHD) simulations of the MRI. 

With the disk's highly conducting surface neglected, the 
fast outward diffusive drift of the large-scale poloidal mag- 
netic field in a turbulent disk can be readily understood 
by looking at the induction equation for the vertical field 
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{rBz)u — Tjr 
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where u > is the accretion speed and rj is the turbulent 
diffusivity both assumed independent of z (Lovelace et al. 
1994). (Both assumptions are invalid and are removed in 
this work.) For a dipole type field, B^ is an even function 
of z and it changes very little over the half-thickness h of 
a thin disk {h ^ r). Consequently, the average of this 
equation over the half-thickness gives 
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where {Br)h is the radial field at the disk's surface. The 
terms inside the square brackets represents the flux of the 
z— magnetic flux. For stationary conditions it vanishes. 
The term proportional to u describes the inward advec- 



tion of the magnetic field, the term oc rj{Br)h describes 
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the outward diffusive drift of the field, and the last term 
represents the radial diffusion of the field. For {Br)h ~ 
the radial difi^usion term is negligible compared with the 
diffusive drift term. For a weak field and the standard 
a— disk model with turbulent viscosity of the order of the 
turbulent diffusivity, the accretion speed is it ~ acgh/r 
(with Cs the sound speed), whereas the outward drift speed 
of the field is ~ acg which is a factor r/h^ 1 larger than 
u. That is, a stationary solution with i?z ^ is not possi- 
ble. A stationary or growing B^ field is possible if the field 
is suflSciently strong to give appreciable outfiow of angular 
momentum to jets (Lovelace et al. 1994). 

Three-dimensional MHD simulations have been per- 
formed which give some information on the advec- 
tion/diffusion of a large scale field. These simulations re- 
solve the largest scales of the MRI turbulence and therefore 
self-consistently include the turbulent viscosity and diffu- 
sivity. Most simulations performed to date have investi- 
gated conditions in which the accreting matter contains no 
net magnetic fltix and where no magnetic field is supplied 
at the boundary of the computational domain (e.g., Hirose 
et al. 2004; De Villiers et al. 2005; Hawley & Krolik 2006; 
McKinney & Narayan 2007) . In these simulations stretch- 
ing of locally poloidal field lines in the initial configuration 
leads to a large-scale poloidal fields and jet structures in 
the inner disk. However, in simulations by Igumenshchev, 
Narayan, & Abramowicz 2003; Igumenshchev 2008) weak 
poloidal flux injected at the outer boundary is clearly ob- 
served to be dragged into the central region of the disk, 
leading to the buildup of a strong poloidal magnetic field 
close to the central object. This fiux build up and its limit 
are discussed by Narayan, Igumenshchev, & Abramowicz 
(2003). The extent to which the magnetic field advection 
seen in numerical simulations depends on having a thick 
disk or nonaxisymmetric conditions is unclear. 

In this work we calculate the profiles through the disk of 
stationary accretion flows (with radial and azimuthal com- 
ponents), and the profiles of a large-scale, weak magnetic 
field taking into account the turbulent viscosity and dif- 
fusivity due to the MRI and the fact that the turbulence 
vanishes at the surface of the disk. By a weak field we 
mean that the magnetic pressure in the middle of the disk 
is less than the thermal pressure. Related calculations of 
the disk structure were done earlier by Konigl (1989), Li 
(1995), Ogilvie and Livio (2001) but without taking into 
account the absence of turbulence at the disk's surface. 
Recent work calls into question the a-description of the 
MRI turbulence in accretion disks and develops a closure 
model which fits shearing box simulation results (Pessah, 
Chan, & Psaltis 2008). Analysis of this more complicated 
model is deferred a future study. 

Section 2 develops the equations for the flow and mag- 
netic field in a viscous diffusive disk. Section 3 discusses 
the boundary conditions at the surface of the disk. Sec- 
tion 4 derives the internal flow/field solutions for an ana- 
lytically soluble disk model. Section 5 discusses the exter- 
nal flow/field solutions which may be magnetocentrifugal 
winds or electromagnetic outflows. Section 6 gives the 
conclusions of this work. 



2. THEORY 



We consider the non-ideal magnetohydrodynamics of a 
thin axisymmetric, viscous, resistive disk threaded by a 
large-scale dipole-symmetry magnetic field B. We use a 
cylindrical (r, (p, z) inertial coordinate system in which the 
time-averaged magnetic field is B = B^-r + B^4> + B^z, 
and the time-averaged flow velocity is v = Vrr+v^<p +VzZ. 
The main equations are 



P-TT = -Vp + pg + -J xB + F"" , 
at c 

9B 

— = Vx (v X B) - V X (r?V X B) . 



(1) 
(2) 



These equations are supplemented by the continuity equa- 
tion, V(pv) = 0, by V X B = 4tt3/c, and by V • B = 0. 
Here, r] is the magnetic diffusivity, F'' = — V-T'^ is the vis- 
cous force with Tj^^, = — pv{dvj / dxk+dvk/ dxj — {2/21)5 ■ 
v) (in Cartesian coordinates), and v is the kinematic vis- 
cosity. For simplicity, in place of an energy equation we 
consider the adiabatic dependence pec p'^ , with 7 the adi- 
abatic index. 

We assume that both the viscosity and the diffusivity 
are due to magneto-rotational (MRI) turbulence in the 
disk so that 

(3) 



v = Vr, = a-f- g{z) , 

where V is the magnetic Prandtl number of the turbulence 
assumed a constant of order unity (Bisnovatyi-Kogan & 
Ruzmaikin 1976), a < 1 is the dimensionless Shakura- 
Sunyaev (1973) parameter, c^o is the midplane isothermal 
sound speed, fix = (GM/r^)^/^ is the Keplerian angu- 
lar velocity of the disk, and M is the mass of the central 
object. The function g{z) accounts for the absence of tur- 
bulence in the surface layer of the disk (Bisnovatyi-Kogan 
& Lovelace 2007; Rothstein & Lovelace 2008). In the body 
of the disk g = 1, whereas at the surface of the disk, at 
say Z5, g tends over a short distance to a very small value 
^ 10~*, effectively zero, which is the ratio of the Spitzer 
diffusivity of the disk's surface layer to the turbulent dif- 
fusivity of the body of the disk. At the disk's surface the 
density is much smaller than its midplane value. 

We consider stationary solutions of equations (1) and (2) 
for a weak large-scale magnetic field. These can be greatly 
simplified for thin disks where the disk half-thickness, of 
the order oi h = Cso/^k, is much less than r. Thus we 
have the small parameter 

h 



£ = 



VK 



<€. 1 



(4) 



It is useful in the following to use the dimensionless height 

C s z/h. 

The three magnetic field components are assumed to 
be of comparable magnitude on the disk's surface, but 
Br = = B^ on the midplane. On the other hand the 
axial magnetic field changes by only a small almount go- 
ing from the midplane to the surface, AB^ ^ sBr <C Bz 
(from V • B = 0) so that B^ « const inside the disk. As 
a consequence, the dBj/dr terms in the magnetic force in 
equation (1) can all be dropped in favor of the dBj/dz 
terms (with j = r, 0). Thus, we neglect the final term 
of the induction equation given in the Introduction. It is 
important to keep in mind that Bj is the large scale field; 
the approximation does not apply to the small-scale field 
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which gives the viscosity and diffusivity. The three velocity 
components are assumed to satisfy v'^ <^ c% and ■^v'^- 
Consequently, (r, z) is close in value to the Keplerian 
value VK{r) = {GM/ry/^. Thus, dv^/dr = -(l/2)(i;^/r) 
to a good approximation. 

With these assumptions, the radial component of equa- 
tion (1) gives 

The dominant viscous force is = —dT^^/dz with 

t;; = -pvdvr/dz. 

We normahze the field components by Bo = Bz{r,z = 
0), with br = Br/Bo, 60 = B^/Bo, and b, = B,/Bo « 1. 
Also, we define = v^{r,z)/vK{r) and the accretion 
speed Ur = — Vr/icxCgo). For the assumed dipole field 
symmetry, br and b^ are odd functions of C whereas Ur 
and U0 are even functions. 

Equation (5) then gives 



dbr 0p 
5C " ~ 



(1 



where p = p(r, z)/ po with po = p{t, z = 0). The midplane 
plasma beta is 



/3 



47rpoc: 



Bl 



(7) 



where fc„ 



-(91n]3/91nr is assumed of order unity 



and p = pCg. Note that /? = c^o/v\q, where vao = 
i?o/(47rpo)^^^ is the midplane Alfven velocity. The rough 
condition for the MRI instability and the associated tur- 
bulence in the disk is /? > 1 (Balbus & Hawley 1998). In 
the following we assume /? > 1, which we refer to as a weak 
magnetic field. 

The 0— component of equation (1) gives 



db. 



(8) 



where ki, = 91n(pCgQr^//i)/51n(r) > is of order unity. 
In addition to the well- know viscous force [F^{a) = 

-r-^d{r'^T^^)/dr with T^^ = -pvrd{v^lr)/dr] which 
gives the term oc fc^, we must include the force contri- 
bution F^{b) = -dT^Jdz with T^^ = -pvdv^jdz. This 
gives the second derivative term in equation (8). 

Note that integration of equation (8) from C = (where 
6^ = and du^jdC, = 0) to Cs + e (where 5 = 0) gives 

6^5+ = ^a/3S(3£A;^ - Ur) , 

where Ur = J^'^ dC, pur/Ti is the average accretion speed, 

E = d( p, and the S+ subscript indicates evaluation 
at C = Cs + The average accretion speed, written as 



Ur = Uo — 



af3T, 



(9) 



is the sum of a viscous contribution, uq = 3ek,y, and a 
magnetic contribution (oc b^s+) due to the loss of angular 
momentum from the surface of the disk where necessarily 
b4>s+ < (Lovelace, Romanova, & Newman 1994). Equa- 
tion (9) is discussed further in §5. The continuity equation 
implies that rhpoT,{aCsoUr) is independent of r. 



The z— component of equation (1) gives 



(10) 



The neglected viscous force in this equation, 
r~^d{rTrz)/dr with Trz = —pvdVr/dz, is smaller than 
the retained pressure gradient term by a factor of order 
a^e <C 1. The term involving the magnetic field describes 
the magnetic compression of the disk because b'^+b'^ at the 
surface of the disk is larger than its midplane value which 
is zero (Wang, Sulkanen, & Lovelace 1990). For /3 » 1 
the compression effect is small and can be neglected. 

As mentioned we assume p (x p^ which can be written 
^P = Pcio{p/poV~^- Thus 

for /3 » 1. The density goes to zero at Cm = [27/(7 ^ 
1)]^/^. However, before this distance is reached the MRI 
turbulence is suppressed, and g{C) in equation (3) is effec- 
tively zero. 

The toroidal component of Ohm's law (equivalent to 
equation 2), = a{v x B)^, with a = (? /{A-kt}), gives 

Multiplying this equation by g, integrating from ^5 — e to 
(s + e, and using the facts that dg/d( = —S{( — (s) and 
that \ur\ is bounded implies that brs+ = Ws-- That is, 
there is no jump in br across the highly conducting surface 
layer. 

Note that multiplying equation (12) by g and integra- 
tion from C = to + e gives 

brS = VCs{Ur) , (13) 

Where (..) = /^^^ dC(-)/C5. 

The other components of Ohm's law give 

duA, 3e ae d f db^\ 

Combining equations (6) and (12) gives 

ur=^-§{l-k,e^-ul) ■ '^'^^^ 

1 -I- 6U4, 



«^ I . (15) 



f 1" ' V dC 
For thin disks, e ^ 1, and /? > 1, we have u,^ 
with {Su^Y <^ 1 which follows from the integral of equa- 
tion (6). Consequently, 

-^"^ = - Si + "wdc [''"w) ' ^''^ 

to a good approximation. 

We first take the derivative of equation (14) and then 
substitute the br derivative with equation (12). In turn, 
the U0 derivatives can be put in terms of Ur and its deriva- 
tives using equation(16). In this way we obtain 

4.2 f d d fi d / dur 



52 



d 



. d 



/id_ 



tlr_ 

pg 

dUr 



+ ^'0'^{p9{ur-9Uo))+P'^l^l 



+ 3pP^^ = . (17) 
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The equation can be integrated from C = out to the 
surface of the disk (s where boundary conditions apply. 
Because Ur is an even function of C, the odd derivatives 
of Ur are zero at ^ = and one needs to specify Ur(0), 
u"(0), and m*"(0). A "shooting method" can be apphed 
where the values of Ur{0), '«"(0), and u™{0) are adjusted 
to satisfy the boundary conditions. Once Ur{C) is calcu- 
lated, equations (8), (12), and (14) can be integrated to 
obtain 60(C), ^^.(C), and u^{0- 
For specificity we take 

> (18) 



where Cs < Cm and 5^1. That is, we neglect the ratio 
of the Spitzer diffusivity on the surface of the disk to its 
value in the central part of the disk. An estimate of Cs can 
be made by noting that /3(C) = 4ttp{()/B^ = /3(p)'^ w 1 at 

Cs. This gives Ci/Cm = 1 -/3-(^-i)/T and ps/po = P'^'^. 

If £ denotes the fraction of the disk mass accretion rate 
which goes into outflows, then we can estimate the ver- 
tical speed of matter at the disk's surface as Wz(Cs) ~ 
£{2h/r){po/ps)\vr\ = E{2h/r)l3^/^\vr\. Our neglect of 
evidently requires that E{2h/r)l3^/^ <^ 1. 

3. BOUNDARY CONDITIONS 

We restrict our attention to physical solutions which (a) 
have net mass accretion, 

M = AirrhpoaCsoT, Ur > , (19) 
and (b) have 6^ < on the disk's surface. This condition 
on b^s+ corresponds to an efflux of angular momentum 
and energy (or their absence) from the disk to its corona 
rather than the reverse. From equation (9), the condition 
bct,s+ ^ is the same as Ur > uo, where uo is the minimum 
(viscous) accretion speed. Note that Ur/uo — 1 is the frac- 
tion of the accretion power which goes into the outflows or 
jets (Lovelace et al. 1994). Clearly, the condition on 6^5+ 
implies that M > so that there is only one condition. 

In general there is a continuum of values of 5: 

for the considered solutions inside the disk. The value of 
b^s+ can be determined by matching the calculated fields 
brs and 6<jis+ onto an external field and flow solution as 
discussed in §5. 

From equation (12) we found that there is no jump in 
br across the conducting surface layer. Thus, integration 
of equation (6) from Cs — e to Cs + e implies that 

This represents a second condition on the disk solutions. 
Integration of equation (14) from Cs ^ e to Cs + e gives 

ae db^ 

This velocity jump must be zero as can be shown by in- 
spection of the total angular momentum flux-density in 
the 2:— direction, rT^z = —Bfj^Bz/Ai: — pvidv^jdz), where 
the first term is the magnetic stress and the second is the 
viscous stress. A jump in would give a delta-function 
contribution to the viscous stress which cannot be bal- 
anced by the magnetic stress. Therefore 

^^'^ = . (21) 
Cs- 



9C 



This is a third condition on the disk solutions. 

Integration of equation (8) from Cs ^ e to Cs + e gives 



'0S+ 



50S- 



dC 



(22) 



This equation is equivalent to the continuity of the an- 
gular momentum flux-density across the surface layer, 
namely, rT^^ = —B'^Bz/'^'n (above the layer) is equal to 

rT-^ = -B-B,/ATT-pv{dv^/dz)\. (below the layer) . The 
jump bfj,s+ — b^s~ corresponds to a radial surface current 
flow in the highly conducting surface layer of the disk, 
Jr = -(c/47r)Bo(605+ - 60S-)- 
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Fig. 1. — Radial and toroidal field components (normalized to B^) 
at the disk's surface as a function of the average accretion speed 
Ur (normalized by the viscous accretion speed uq). For this plot 
(3 = 100 and Prandtl numbers V = 1 and 2. Note that is 
given by equation (10) and is independent of V and b^s is given by 
equation (13). 



4. INTERNAL SOLUTIONS 

Here, to simplify the analysis we consider the limit 
where 7 — s- 00 in equation (11) and ^ — > in equation (18). 
Then, Cs ~^ Cm and both p and g are unit step functions 
going to zero at Ct?j. = V2. Also, Ur — (ur) and S = \/2. 
Thus the above physical condition Ur > Sefe^ = uq implies 
that brs > uqCs'P from equation (13). We assume kp and 
K = 1. 

The solutions to equation (17) are Ur oc exp(i/cjC) (with 
j — 1, 2, 3), where 

a'^(5^{k]f -t- 2Va^(3{k]f + {a^ (3^ + V^)k] ~ SPV"^ = , 

(23) 

is a cubic in fc| . The discriminant of the cubic is negative 
so that there is one real root, kf, and a complex conjugate 
pair of roots, 3. Because Ur is an even function of C we 
can write 



— ai cos(fcoC) + ^2 cos(fcrC) cosh(fciC) 
-I- 03 sin(fcrC) sinh(fciC) , 



(24) 



where ko ~ \/f^, kr = Re(-\/fcf), and ki — Im(-\/A!f). 

The three unknown constants ai, 02,0,3 are reduced to 
two by imposing equation (20). The two are then reduced 
to one constant by imposing equation (21). The remaining 
constant is restricted by the condition b^s+ 5: 0. 

We consider a thin disk, e = h/r = 0.05, and a viscosity 
parameter a = 0.1. Figure 1 shows the dependences of the 
surface field components on the average accretions speed 
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for P = 100 and two values of V. The b^s dependence is 
given by equation (9) and is independent of P while the 
brs is given by equation (13). 




(r-ro)/h 

Fig. 2. — Radial flow speed Vr = —Mr (normalized to aCgo) as a 
function of = z/h and a sample poloidal (Br , Bz ) magnetic field 
line for /3 = 10^ and "P = 1. Also, ro is a reference radius. 



surface layer ^ 




Fig. 3. — Toroidal magnetic field = B^/B^ and toroidal veloc- 
ity <5m0 = (d0 — vk)/vk (with vk the Keplerian velocity) for the 
case where /3 = 100 and V = 1. The jump in the toroidal magnetic 
field at the disk's surface is shown by the dashed line. 

Figure 4 shows the profile of the accretion speed Mr(C) 
and the shape of the poloidal magnetic field Bp for (3 = 300 
and V = \. Figure 5 shows the profiles of the toroidal mag- 
netic field hff, and the fractional deviation of the toroidal 
velocity from the Keplerian value, 5uip = {v^ — vk)/vk- 
For this case, Ur/uQ — 1.5, b^s+ — —1.59, and brS = 0.318. 
Inspection of the flow/field solution shows that the angu- 
lar momentum flux rT^z > in the top half of the disk 
and at the disk's surface this flux goes into an outflow or 
jet. 



Figure 2 shows both the profile of the accretion speed 
Ur(C) and the shape of the poloidal magnetic field Bp 
for /3 = 100 and V — \. Figure 3 shows the profiles 
of the toroidal magnetic field b^ and the fractional devi- 
ation of the toroidal velocity from the Keplerian value, 
5U(f, — {vcj, — vk)/vk- For this case, Ur/uo — 1.30, 
b^s+ = —0.321, and brs = 0.276. In this case (and a 
range of others discussed below) , we find that there is ra- 
dial inflow of the top and bottom parts of the disk whereas 
there is radial outflow < of the part of the disk around 
the midplane. Inspection of the flow/fleld solution shows 
that the top and bottom parts of the disk lose angular 
momentum (by the vertical angular momentum flux rT^z ) 
both to (a) magnetic winds or jets from the disk's surfaces 
and to (b) the vertical flow of angular momentum to the 
midplane part of the disk which flows radially outward. 

We find that the flow pattern is the same as in Figure 
2 for 10 < /3 < 200 and V > 1. For larger values of /3 and 
V — 1, the flow pattern changes from that in Figure 2 to 
that in Figure 4 for f3 — 300. However, for /3 = 300 and 
■p = 2, the flow pattern is again similar to that in Figure 
2. For f3 = 10'^ and V = 3 the flow pattern is also the same 
as in Figure 2. For smaller viscosity values, a < 0.03, and 
P = 100 there are multiple channels with for example three 
layers of radial inflow and two layers of radial outflow in 
the disk. 
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Fig. 4. — Radial flow speed Vr = —Ur (normalized to ocso) as a 
function oi C, = z/h and a sample poloidal (Br,Bz) magnetic field 
line for /3 = 300 and P = 1. Also, ro is a reference radius. 



5. EXTERNAL SOLUTIONS 

As mentioned in §3, the value of 60s+ < is not fixed 
by the solution for the field and flow inside the disk. Its 
value can be determined by matching the calculated sur- 
face fields brS and t^s-i- onto an external magnetic wind or 
jet solution. Stability of the wind or jet solution to current 
driven kinking is predicted to limit the ratio of the toroidal 
to axial magnetic field components at the disk's surface 
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\b^S+\ to values ^ 0{2Trr/L^) (Hsu & Bellan 2002; Naka- 
mura, Li, & Li 2007), where is the length-scale of field 
divergence of the wind or jet at the disk surface. From 
known wind and jet solutions we estimate 2TTr/Lz ~ tt 
(Lovelace, Berk, & Contopoulos 1991; Ustyugova et al. 
1999; Ustyugova et al. 2000; Lovelace et al. 2002). Re- 
call that Ur/uo — 1 = 2\bij,s+\/ (aPJ^uo) (from equation 9) 
is the faction of the accretion power going into the jets 
or winds. For the mentioned upper limit on 16^5+ 1, we 
find Ur/uQ — 1 <^ C'[27r/(a/3Euo)]. From equation (13) we 
have brs = {'P(suo){{ur) /uq) . Therefore, for /3 ^ 1 and 
(ur) « uo, we have brS ~ VCsuq- 

The matching of internal and external field/flow solu- 
tions has been carried out by Konigl (1989) and Li (1995) 
for the case of self-similar [^^(r, 0) ~ r~^/^] magnetocen- 
trifugally outflows from the disk's surface. These outflows 
occur under conditions where the poloidal field lines at the 
disk's surface are tipped relative to the rotation axis by 
more than 30° which corresponds to brs > 3~^/^ w 0.577 
(Blandford & Payne 1982). The outfiows typically carry 
a significant mass flux. For the internal ficld/fiow solu- 
tions discussed in §4 with f3 ^ 1, we conclude that brs is 
sufficiently large for magnetocentrifugal outflows only for 
turbulent magnetic Prandtl numbers, V ^ 2.7. Shu and 
collaborators (e.g., Cai et al. 2008, and references therein) 
have developed detailed 'X-wind' models which depend on 
the disk having Prandtl numbers larger than unity. Re- 
cent MHD simulations by Romanova et al. (2008) provides 
evidence of conical or X-wind type outflows for Prandtl 
numbers > 1. 
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Fig. 5. — Toroidal magnetic field = B^/ and toroidal veloc- 
ity Su^ = {v^ — vx)/vk (with vx the Keplerian velocity) for the 
case where f} = 300 and "P = 1. The jump in the toroidal magnetic 
field at the disk's surface is shown by the dashed line. 

For Prandtl numbers say V < 2.7, the values of brs are 
too small for there to be a magnetocentrifugal outflow. In 
this case there is an outflow of electromagnetic energy and 
angular momentum from the disk (with little mass out- 
flow) in the form of a magnetically dominated or 'Poynting 
flux jet' (Lovelace, Wang, & Sulkanen 1987; Lovelace, et al. 
2002) also referred to as a 'magnetic tower jet' (Lynden- 
Befl 1996, 2003). MHD simulations have established the 
occurrence of Poynting-flux jets under different conditions 
(Ustyugova et al. 2000, 2006; Kato, Kudoh, & Shibata 
2002; Kato 2007). Laboratory experiments have allowed 
the generation of magnetically dominated jets (Hsu & Bel- 
lan 2002; Lebedev et al. 2005). 



6. DISCUSSION 

A study is made of stationary axisymmetric accretion 
flows [vr{z),V(j,{z),Vz — 0] and the large-scale, weak mag- 
netic field [Br{z), B^(z), Bz ~ const] taking into account 
the turbulent viscosity and diffusivity due to the MRI 
and the fact that the turbulence vanishes at the surface 
of the disk as discussed by Bisnovayi-Kogan & Lovelace 
(2007) and Rothestein & Lovelace (2008). We derive a 
sixth-order differential equation for the radial flow veloc- 
ity Vr{z) (equation 17) which depends mainly on the ra- 
tio of the midplane thermal to magnetic pressures > 1 
and the Prandtl number of the turbulence V = viscos- 
ity /diffusivity. Boundary conditions at the disk's surfaces 
take into account the outflow of angular momentum to 
magnetic winds or jets and allow for current flow in the 
highly conducting surface layers. In general we find that 
there is a radial surface current but no toroidal surface cur- 
rent. The stability of this surface current layer is unknown 
and remains to be studied. If the layer is unstable to kink- 
ing this may cause a thickening of the current layer. We 
argue that stability of the winds or jets will limit the ratio 
of the toroidal to axial field at the disk's surface |60s+| to 
values <^ TT. The stationary solutions we find indicate that 
a weak {(i ^ 1), large-scale field does not diffuse away as 
suggested by earlier work (e.g., Lubow et al. 1994) which 
assumed brs > 3^^^^. 

For a wide range of parameters /3 > 1 and V > I, we 
find stationary channel type flows where the flow is radi- 
ally outward near the midplane of the disk and radially 
inward in the top and bottom parts of the disk. Solutions 
with inward flow near the midplane and outflow in the top 
and bottom parts of the disk are also found. The solutions 
with radial outflow near the midplane are of interest for the 
outward transport of chondrules in protostellar disks from 
distances close to the star (^ 0.05 AU) (where they are 
melted and bombarded by high energy particles) to larger 
distances (> 1 AU) where they are observed in the Solar 
system. Outward transport of chondrules from distances 
~ 0.05 AU to > 1 AU by an X-wind has been discussed 
by Shu et al. (2001). 

The flow/fleld solutions found here in a viscous/diffusive 
disk and are different from the exponentially growing chan- 
nel flow solutions found by Goodman and Xu (1994) for 
an MRI in an ideal MHD unstable shearing box. Chan- 
nel solutions in viscous/diffusive disks were found earlier 
by Ogilvie & Livio (2001) and by Salmeron, Konigl, & 
Wardle (2007) for conditions different from those consid- 
ered here. In general we find that the magnitude of the 
toroidal magnetic field component inside the disk is much 
larger than the other field components. The fact that the 
viscous accretion speed is very small, aecso, nieans that 
even a small large-scale field can significantly influence the 
accretion flow. We flnd that Prandtl numbers larger than 
a critical value estimated to be 2.7 are needed in order for 
there to be magnetocentrifugal outflows from the disk's 
surface. For smaller V, electromagnetic outflows are pre- 
dicted. Owing to the stability condition, |6^s+| ^ tt, the 
fraction of the accretion power going into magnetic out- 
flows or jets is <^ const /3~^ ~ 

Analysis of the time-dependent accretion of the large- 
scale B— fleld is clearly needed to study the ampliflcation 
of the fleld and build up of magnetic flux in the inner region 
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of the disk. One method is to use global 3D MHD sim- 
ulations (Igumenshchev et al. 2003; Igumenshchev 2008), 
but this has the difhculty of resolving the very thin highly 
conducting surface layers of the disk. Another method is 
to generalize the approach of Lovelace ct al. (1994) taking 
into account the results of the present work. This is pos- 
sible because the radial accretion time {r/\ur\) is typically 
much longer than the viscous diffusion time across the disk 
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